Glassy correlations and microstructures in randomly crosslinked homopolymer blends 
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We consider a microscopic model of a polymer blend that is prone to phase separation. Permanent 
crosslinks are introduced between randomly chosen pairs of monomers, drawn from the Deam- 
Edwards distribution. Thereby, not only density but also concentration fluctuations of the melt 
are quenched-in in the gel state, which emerges upon sufficient crosslinking. We derive a Landau 
expansion in terms of the order parameters for gelation and phase separation, and analyze it on 
the mean-field level, including Gaussian fluctuations. The mixed gel is characterized by thermal as 
well as time-persistent (glassy) concentration fluctuations. Whereas the former are independent of 
the preparation state, the latter reflect the concentration fluctuations at the instant of crosslinking, 
provided the mesh size is smaller than the correlation length of phase separation. The mixed gel 
becomes unstable to microphase separation upon lowering the temperature in the gel phase. Whereas 
the length scale of microphase separation is given by the mesh size, at least close to the transition, 
the emergent microstructure depends on the composition and compressibility of the melt. Ifexagonal 
structures, as well as lamellae or random structures with a unique wavelength, can be energetically 
favorable. 

PACS numbers: 82.70.Gg, 64.75.-l-g, 61.43.-j 
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I. INTRODUCTION 

Crosslinked homopolymer blends exhibit a rich phase 
diagram, due to the competition between phase separa- 
tion and crosslinking. The simplest case is of a blend of 
two homopolymer species, "j4" and "B", whose incompati- 
bility is controlled by the Flory-Huggins parameter x a-nd 
which are crosslinked irreversibly by some number of 
chemical bonds. In addition, the concentration fluctua- 
tions can be controlled independently in the process of 
crosslinking (preparation state) and the well crosslinked 
gel (measurement state), e.g. by lowering the tempera- 
ture in the gel. Hence we have three control parameters: 
the incompatibility Xp in the preparation state, the in- 
compatibility Xm in the measurement state, and the num- 
ber of crosslinks per chain, fi — Nc/N, where N denotes 
the total number of chains in the melt. 

A statistical mechanical theory thus has to include not 
only the average over the quenched disorder (crosslink re- 
alization) but also the "memory" of the preparation state. 
This can be achieved in the following way: We start 
from a microscopic model, which accounts for the repul- 
sive interaction of all monomers, irrespective of species 
(excluded volume), as well as for a repulsive interac- 
tion between the different species only (incompatibility) . 
Crosslinks are introduced between randomly-chosen pairs 
of monomers. The probability for a particular crosslink 
configuration depends on the preparation state of the sys- 
tem, such that monomers with a high probability to be 
close in the preparation state have a high probability to 
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be crosslinked. Thereby the crosslinks indeed preserve 
the memory of the preparation state. Mathematically, 
this is achieved via the Deam-Edwards distribution Q| 
and the replica trick to average over the quenched disor- 
der. 

We expect and indeed find signatures of the prepara- 
tion state in the gel. An example are the concentration 
fluctuations which are frozen in by the crosslinks. If the 
preparation state is close to macroscopic phase separation 
then the glassy, i.e. time persistent, concentration fluctu- 
ations reflect the correlations of the melt at the moment 
of crosslinking. If, on the other hand, the preparation 
state is far from phase separation then the frozen fluctu- 
ations are completely random and follow the pattern set 
by the crosslinks. 

Lowering the temperature in the gel or, equivalently, 
increasing the incompatibility at measurement Xm will ul- 
timately give rise to microphase separation, while macro- 
scopic demixing is suppressed by the crosslinks. A vari- 
ety of microphases can exist, depending on the compo- 
sition of the blend and its compressibility. If the mix- 
ture is symmetric, having equal concentrations of and 
"5" monomers, then lamellae are energetically favorable, 
whereas for an asymmetric mixture hexagons prevail. 
A finite compressibility enhances the tendency towards 
phase separation and can induce a random pattern, con- 
sisting of a superposition of many lamellar phases of dif- 
ferent orientation. In all cases the critical wave-number 
is given by the mesh size or localization length of the gel. 

The subject of crosslinked homopolymer blends was 
first addressed by de Gennes who pointed out that 
the mixed state is stabilized in the gel and eventually 
undergoes microphase separation. His predictions were 
verified experimentally ^ with, however, a discrep- 
ancy in the scattering intensity for small wave-number. 
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This was traced back to the neglect of concentration 
fluctuations which are present during crosslinking and 
are partially frozen in by the crosslinks. Subsequently 
several attempts were made to include these effects ap- 
proximately 1^, H, H, 1^]. Studies of crosshnked systems, 
based on the microscopic model by Panyukov and Ra- 
bin were reported by Sfatos and Shakhnovich Q; as 
far as homopolymer blends are concerned, these authors 
recover de Gennes' result within a microscopic approach. 
Computer simulations were carried out by Lay and Som- 
mer who studied, in particular, the relation of the 
domain sizes to the mesh size of the gel. In the last sec- 
tion of our paper we present a detailed discussion of the 
literature in comparison with our own results. 

The paper is organized as follows: In sec. ^ we formu- 
late a microscopic model of crosslinked polymer chains. 
Subsequently (sec. Ill) we derive a Landau expansion in 



terms of the order parameters for gelation and phase sep- 
aration. The Landau theory allows us to discuss the 
mixed gel (sec. IV) as well as microphase separation 
(sec. 0). We conclude with a short summary, a com- 
parison with previous theoretical work and an outlook. 
A short account of our results for the special case of an 
incompressible melt with equal concentrations of A and 
B monomers was given previously in ref . . 



all lengths are measured in units of the radius of gyra- 
tion of the free chains, Rg = Lb^ /6. The rescaled volume 
reads V := {d/Lb^Y^^ V. 

The chain connectivity is described by the usual 
Wiener Hamiltonian 



H 



Na 



a=A,B i=l 



ds 



(1) 



the excluded volume term controlling compressibility 
reads 



xS{ra4s)-ra'^^'is')), (2) 



and the incompatibility of the two monomer species is 
modeled by the interaction 



VX 



Na ^1 ^1 

i^Sa.a' - 1) ^ ds ds' 

a,a'=A,B i,i' = l "'0 JQ 



^S{raAs)-ra'^Hs')) . (3) 



II. MODEL 

A. Uncrosslinked homopolymer blend 

We first consider an uncrosslinked blend of polymer, 
modeled as a system of Gaussian phantom chains of equal 
degree of polymerization L and step length b. The melt 
is taken to contain Na chains of type A and Nb chains 
of type B. In general, there will be an imbalance in con- 
centration q {Na - Nb)/N where N ^ Na+ Nb \s 
the total number of chains, occupying a volume V in d- 
dimensional space. The monomers positions are denoted 
as Ra,i{s), where a — A,B refers to the chain species, 
i — 1, . . . , Na enumerates the chains, and s = ... 1 is 
the continuous index for a "site" on a chain. It turns out 
to be convenient to express the monomer positions by 
dimensionless vectors T"a,i(s) = ^fdfUP- RaA{s), so that 



Although the chain elasticity and the volume exclusion 
are of mainly entropic origin, the incompatibility is as- 
sumed to be a chiefly energetic contribution. Neverthe- 
less we let fceT = 1 in the following to simplify the ex- 
pressions. Instead of changing the temperature, we shall 
tune A and particularly x, which will serve as the inverse 
temperature. 



B. Crosslinking 

Chemical crosslinking induces a random number M 
of permanent bonds between randomly selected pairs of 
monomers; a particular realization of crosslinks is de- 
noted by C = {(fle, le, SejOg, Zg, Sg)}g£]^. The links are 
modeled as hard constraints with zero bond length. The 
partition function of the crosslinked melt, relative to a 
melt of non-interacting chains, thus reads 



M 



W 



Z{C) := / []^(r„^,,,(sg)-ra,,,,(s'g))exp{-ff^-ijx} 

\ e=l 

_j'Dra.,{s) neli^K,»e(ge)-r,>^,,(4)) CXp {-g^ - - } 

/2?r,,,(s) exp{-i/^} 

I 



(4) 



Here we have implicitly defined the expectation value crosslinked melt. 
(. . . )^ with respect to the Hamilton function of the un- 
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C. Disorder average and Deam-Edwards 
distribution 

We are interested in the properties of the "generic" melt 
rather than in the properties of a melt with a specific set 
of crosslinks. Furthermore, we assume the system to be 
self-averaging in the thermodynamic limit. Therefore we 
will consider disorder averages of the observables with 
respect to the quenched randomness of crosslinks. 

We specify the probability distribution of the crosslink 
sets following the strategy of Deam and Edwards |l| . We 
suppose that the dominant crosslink sets are those that 
are most compatible with the uncrosslinked melt. More 
precisely, we assume a probability distribution 

Pm(C)(x Zp(C). (5) 

Here, is given by eq. (Q), evaluated at A = Ap and x = 
Xp, which characterize the system prior to crosslinking. 
Disorder averages with respect to W'm will be denoted by 
square brackets. 

D. Order parameters for the homopolymer blend 

To discriminate between the liquid state and the amor- 
phous solid state of the polymer system we use the order 
parameter proposed in [l^ : 

f^/ei ... fc, := ^ / ds(exp(ifcira,i(s)))^ 

a=A,B 1=1 •'^ 

X • • • X (exp (ifcgra,i(s)))^ , (6) 

for g = 1,2,... and nonzero {k^}. The symbol {■ ■ ■)c 
denotes the thermal expectation value in the presence of 
a particular realization C of crosslinks. In the case 5 = 1, 
eq. (H) is the thermal average of the monomer density in 
Fourier space, 

^^-^ E E A.e---(^). (7) 

In the liquid state, a monomer explores the sam- 
ple volume uniformly. Hence, the equilibrium value of 
the local density is constant and the Fourier transform 
{exp{ikri{s)))c vanishes (except for A; = 0, which we ex- 
clude). The order parameter therefore is always zero 
in the liquid state. 

In a solid, at least a finite fraction of the monomers 
are localized about points ba,i{s) in space. For these 
monomers, (exp(iA;ra,i(s)))c oc e''^''^-'^'') ^ 0. However 
for an amorphous, i.e. macroscopically translationally 
invariant (MTI) solid, the disorder averaged expectation 
value [^ki ... kg] vanishes unless ki + - ■ ■ + kg = 0, see Q. 



In particular [{pk)c] = in the MTI state. Hence we can 
discriminate between the liquid and the amorphous solid 
state by means of the g > 2 values of eq. (H). (For the 
signature of crystalline and globular states, see 

Throughout this article, we will also refer to the two 
monomer species as "opposite charges". The identifica- 
tion of ^s and Bs with positive and negative charges, 
respectively, leads to a natural choice for an order pa- 
rameter detecting phase separation: the "charge density" 

^ Na ^1 , Nb „i 

^=^E / d5e-'^-'(^)-lg / dse-'^-.'(^) 

i=l 1=1 

(8) 

measuring the local imbalance of the concentrations of A 
and B. 

In the general case of an asymmetric blend, in which 
there is an excess of either A- or i3-chains, the average 
charge density is given by q ■ N/V — {Na ~ Nb)/V, so 
that ^'fc serves as an order parameter only for fc ^ 0. 
Homogeneous phase separation is indicated by a nonzero 
expectation value of the order parameter in the limit 
k 0. Microstructures, e.g. lamellae or hexagonally or- 
dered cylinders, give rise to a nonzero expectation value 
of [(^'fc)c] at finite wave-number. A nonuniform charge 
density in general is accompanied by mass density mod- 
ulations, except for the incompressible case. 

In the gel phase we expect to find static charge fluc- 
tuations (^fc)c 7^ for all wave-numbers. If the gel 
state is a homogeneous mixture of A and B chains 
then the disorder averaged charge density vanishes, i.e. 
[{^k)c] — 0. The frozen-in fluctuations can only be 
detected by the glassy correlations [(^fc)c(^-fe)c]- In 
general, the quadratic expectation values [{^k^-k)c] 
and [(5'fc)c(^-fc)c] measure volatile and time-persistent 
charge fluctuations in an a priori homogeneous mixture. 

III. EFFECTIVE FREE ENERGY 

The disorder averaged free energy 

\(zV] - 1 

F = -\\nZ(C)]=-lim^^-!-J (9) 

n^O n 

is computed with help of the replica trick. The n-th 
power of Z is made explicit using n independent copies 
- replicas - of the system and an additional replica is 
introduced to account for the Deam-Edwards distribu- 
tion Jl^ . The average over the disorder can be carried 
out explicitly, yet at the cost of a coupling between the 
formerly independent replicas, yielding 

[Z"] =: Z„/Zo , (10) 

where the replicated partition function is given by 
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47V 
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E E 



P n 



0) 



TV 



ri+1 



(11) 



Here, the ^{s) denote the monomer positions in the 
^th j-gpiica,, is the rephcated Wiener average, 

and -ff^+i and H^^^ denote the rephcated Hamiltonians 
of the excluded volume and incompatibility interactions. 
The denominator and the zeroth replica in the numer- 
ator are due to the Deam-Edwards distribution (||) and 
reflect the situation prior to crosslinking. Thus, we have 
to distinguish between the zeroth replica, characterized 
by Ap and Xp (preparation ensemble), and the other n 
replicas, reflecting the situation after crosslinking, char- 
acterized by Am and Xm (measurement ensemble). To 
account for the particular role of the zeroth replica we 
use the notation 



A" := 



Ap if a = 0, 

Am otherwise, 

Xp if a = 0, 

Xm otherwise. 



and (12) 
(13) 



The many-particle problem of the polymer melt can 
be formally reduced to a two-chain problem (one chain 
of each species). It is convenient to introduce {n + 1)- 
fold replicated vectors x :— (a;°, . . . , a;"), and to express 
the exponents in eq. (|l^) in Fourier space, which leads 
to terms quadratic in the monomer and charge densi- 
ties. These can be linearized by means of Hubbard- 
Stratonovich transformations, yielding 



B„- Jv{{^,^,p})exp{~nJ^,,{{^,n,p})} (14) 



with a constant Bn 
l)F-"-l)/i)} -1 



= exp { ^(-"-Am + nq'^Xm + {{n + 
0{n) and the effective free energy 



nM{^,^p}) = f E' ((^ - y 1*^1' + + ^{pt^-. + p-un. 

+ ^EI"fcl'-^^i^'^+(i*'^'/'})-^si^^-(i*'^''^})' (15) 

where 

z±({*,n,p}):-/exp(^'( + ip^-(qTl)*^) Tds e"^'-"^^) + F'" ^% • f As<^^'^A\ . (16) 

I 



The Wiener average now runs over a single replicated 
chain having monomer positions f{s). The auxiliary 
fields \l/fc, pk and 17^, are pairwise dependent via fi ^, — 
(ri j,)* etc., thus the and Vl integrations in ( p^ are 

restricted to the half spaces fc • n > and k ■ h > (with 
arbitrary nonzero constants n and n) . The sums over the 
fe" and k are split into replica sectors: The constant part 
Bn is composed of the fe" = and k = contributions 
(zero replica sector). The symbols '}2'k ^^^"^ X]fc denote 
the sums over nonzero k {single replica sector) and over 
k with nonzero k" in at least two replicas [higher replica 
sector), respectively; both restricted to the above sub- 



spaces. 

As expected, crosslinks give rise to an attractive in- 
teraction between the chains, which, in the absence of 
excluded volume, would cause the chains to collapse into 
a globular state. A sufficiently strong excluded volume 
interaction prevents this collapse, as can be read off the 
coefficient A" :— A° — p/V" of the term quadratic in the 
density. Stability requires Xp > p and Am > p. 

The physical observables, moments of the local density 
and charge density, are related to expectation values of 
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the fields according to 
for fc 7^ and a > 1, and 



and 



for fc = (O, fci, 



i lim(l7,-)-^^, 



,0) with g>2. 
In asymmetric blends, which have an excess of ei- 
ther ^ or B chains, the average charge density qN/V — 
{Na — Nb)/V is nonzero. To simplify the Landau expan- 
sion of the free energy it is then advantageous to either 
work with the fluctuations S'i' = — g of the charge 
density around its mean value or, alternatively, shift the 
monomer density, as is done here. 



method, i.e. by the value of the integrand at the point 
H) making the integrand stationary: 



const • e 



.(*,p,n) 



(17) 



where, by definition, fi) satisfy the stationarity 

conditions 



= , — - 



= 0, 



= 0. 



(18) 



A. Homogeneously mixed liquid state 



IV. HOMOGENEOUSLY MIXED STATES 

On the mean field level, we approximate the functional 
integral over p and 5- in (|l|) by using the saddle point 



One solution of the stationarity conditions ( [l8[ ) is the 
trivial saddle point ^ = p = = Q, corresponding to 
the homogeneously mixed liquid state. To assess its sta- 
bility we consider the Landau expansion to leading order 
around this point: 



2n 
N 



.F„ p, 17}) = ^' (^i, -f 5D (fc^)) I I ' + 5] (i - 5D (P)) I 



a,fc 



2ig^' ^ + E'(^ - - (1 - 9')ffD(fe')) 1*^1'+ o(*^p^ 1^2) , (19) 



r 



where the Debye function 170 is defined in appendix B 1 



The stability limits of the homogeneous liquid can be 
read off from the quadratic coefficients. As g^ik^) de- 
creases monotonically from one to zero, stability against 
solidification and demixing require fj, < 1 and (xp, Xm) < 
1/(1 — g^), respectively. Throughout this article, we as- 
sume that (Ap,Am) > i.e. that the excluded volume 
interaction is strong enough to prevent density instabili- 
ties (see above). 



It should be noted that the condition (Xp,Xm) < 
1/(1 — q^) denote local stability limits only. In mean- 
field theory, the phase separation transition for symmet- 
ric blends is of second order, so the phase transition co- 
incides with the limit of local stability. In asymmetric 
blends the transition is of first order. The loss of local 
stability, as given by the conditions above, then defines 
a spinodal, and the transition occurs at a lower value 
of X- The location of the spinodal depends on the average 
charge q, with a larger critical incompatibility (i.e. lower 
critical temperature) for more asymmetric mixtures. 



The gelation of the homogeneous liquid, driven by in- 
creasing the crosslink concentration, and the microphase 
separation of the resulting gel, induced by cooling, will 
be addressed in the following sections. 



B. Crosslinking in the homogeneously mixed state 

In the liquid state, the polymer blend phase separates 
macroscopically beyond the demixing threshold. The 
subsequent gelation of such a macrophase-separated melt 
would result, apart from the interface, in just two pieces 
of gel having different compositions. It is more interest- 
ing to consider a gel prepared from a homogeneous melt 
to study phase separation in the gel phase. As we shall 
see below, such a gel shows glassy charge density pat- 
terns and, as anticipated, microphase separation instead 
of macroscopic demixing. Therefore, the discussion will 
be restricted to crosslinking in a homogeneously mixed 



blend, Xp < 1/(1 
for q^O. 



q ), including undercooled mixtures 
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Upon gelation, the saddle point = will become 
unstable, making it necessary to complement the expan- 



sion (|l|) of the free energy with the third-order terms: 



f^„({*,p,fl}) = > ' 1 



(20) 



fel,: 



The vertex functions of the cubic terms have been ap- 
proximated by their zero wave-number va lues , the com- 
plete expressions being given in appendix B 1 . This ap- 



proximation is well justified, because the gelation tran- 
sition is always continuous, so that the relevant length- 
scales are very large, compared with the scales of the 
microscopic correlations. Here, we have taken the limit 
of an incompressible melt, which is achieved by integrat- 
ing out the density fluctuations on the Gaussian level and 
subsequently taking the limit X —^ oo. 

We first discuss a gel in the homogeneously mixed state 
= 0, assuming Xm < 1/(1 ^ 9^) • Following Q, we 
consider the order parameter hypothesis 



Q / dTp(T)exp(-fcV2T) 



(21) 



with the shorthand k := X]q=o'''"- Hsre, Q denotes the 
fraction of chains that are localized, i.e. the gel fraction; 
the localization lengths arc distributed according to the 
distribution function p(r). Both have to be determined 
self-consistently as a solution of the stationarity condi- 
tions dH). 

The first two of the stationarity conditions are satisfied 
for any Q and p{t). The third condition is independent 
of the incompatibility parameter. Hence, in the homo- 
geneously mixed regime, the task of determining Q and 
p{t) on the saddle-point level is exactly the same as in 
the pure gelation problem in fl^ . In the present nota- 
tion, the result for the solid state, i.e. ^ > 1, reads 



<5^o./i(M-l)-u; fc72(/i-l) , (22) 



with the gel fraction approximately given by Q ~ 2(/i — 
1). In the liquid state, 0.^—0 and Q — 0. The scal- 
ing function u^x) is defined in (see appendix |C l| 
for details); for convenience we define the shorthand 

w{x) := uj 



C. Stability of the homogeneously mixed gel 

Starting from a gel prepared from a homogeneous melt, 
i.e. Xp < 1/(1 ~ 9^)1 we now allow the incompatibility 
to be changed after crosslinking. In order to keep the gel 
homogeneous, Xm must remain smaller than the critical 
value Xcrit for (micro-)phase separation. As can be seen 
from the term coupling '^k and flf, in the effective free 
energy (|2^), the gel network stabilizes the mixed state; 
the details are discussed in the following. 

As long as the gel is homogeneous, the order parame- 
ter ( p2| ) solves the stationarity conditions (18). To deter- 
mine the stability of the mixed state, we need the second 
derivative of Tn with respect to the charge density, eval- 
uated at the saddle point. We restrict the discussion to 
a weak gel, i.e. /it — 1 <C 1, so that the saddle point value 
of is small and the Hessian can be approximated by 
its expansion to linear order in VU. It then can be read 
off from the Landau expansion (|20|) with fi^ replaced by 
the explicit saddle point value ( |22[ ). The latter vanishes 
unless A; = 0: hence there is no coupling between the 
different fc's, and the Hessian can be calculated indepen- 
dently for each wave-vector. We obtain 



where 



N{l~q^)-A^,^,ik)., 



( c -b ■■■ -b\ 
-ha ■• : 

\-b ... -b a J 



(23) 



(24) 



with 



b■.^^l{fl-l)w{kyifl~l)) 



and 



(25) 



The stability of the homogeneous state is equivalent to 
the positivity of A . In the limit n ^ 0, its eigenvalues 
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FIG. 1: Stability parameter \2{k) for Xm = 1/(1 — q^) as a 
function of k^. 



are given by 

Ai(fc) :=c 
A2(fc) ■.= a + b 



(non-degenerate) , 
(rt-fold degenerate). 



and 



(26) 
(27) 



As we assume crosslinking in the mixed phase, Ai is al- 
ways positive, and thus the stability condition reduces to 
A2 > or, equivalently Xm < Xcrit(M) with 

(1 - 9^)Xcrit(M) 

l/m^x{go{k^) - fi{fi-l)w{kyi^i-l))} . (28) 

Figure |l| shows A2(fc) for Xm = 1/(1 ~ 9^) and different 
crosslink concentrations. Increasing Xm shifts the curve 
downwards. Apparently, an instability towards demixing 
first occurs for a nonzero wave- number kc, which maxi- 
mizes the above expression. 

We consider the case of a weak gel, so /i(/i— 1) ~ /j,— 1, 
and we can write 50 (fc^) ~ l — as the Debye function 
decays much more slowly than w(A;^/(/i — 1)). In this 
approximation, kc is given by 







0X2 
dk^ 



k=k^ 



^w'{kl/{^i-\)), 



(29) 



which leads to k^ ~ 1.61(/i — 1) and (1 — (?^)xcrit — 1 ~ 
kl/Z + (/i - l)wlkl/{^i -1)) « 0.98(m - 1). 

A more precise numerical analysis without these ap- 
proximations yields 



(1 - 9^)Xcrit - 1 = 

0.98 • (Ai - 1) + 0.70 -{ji-lf + {in ~ if) 

and 

fc2 = 1.61-(M-l) + 1.75-(Ai-l)^ + 0((A*- 1)^) 



(30) 



(31) 



The instability for nonzero k implies that the gel un- 
dergoes microscopic (rather than macroscopic) phase sep- 
aration. This is to be expected, because crosslinks per- 
manently connect different chains and thus prevent true 



>< 2 - 








FIG. 2: Stability limits Xcrit(/i) as a function of the asymme- 
try q for the liquid (/i < 1) and gels with different strengths 
{jj, = 1.1 and p — 1.5). 



macroscopic phase separation. The "next best" state for 
the system is phase separation up to the length-scale of 
the network, i.e. the typical mesh size, as given by the 
average localization length f ^ (/i — 1). Hence the insta- 
bility occurs at a critical wave-number kc ~ 1/f. 

The instability is hampered by an increased density fi 
of crosslinks and the asymmetry q of the composition. 
The spinodals for the liquid blend and two solid gels with 
two different degrees of crosslinking are shown in fig. ^ 
as a function of q. In contrast, the critical wave-number 
remains unchanged in agreement with the above argu- 
ment - the critical wave-number is determined by the 
mesh size, which is unaffected by q. The microphase tran- 
sition is addressed in section where it will be shown 
that the average charge also influences the observed mi- 
crostructure. 



D. Pseudo phase diagram 

Three parameters determine the state of the system: fj, 
controls the number of crosslinks, Xp specifies the charge 
fluctuations at preparation and Xm the charge fluctua- 
tions after the gel has been prepared. Each of them 
can be chosen such that the system is close to a criti- 
cal point: fJ. = 1 corresponds to the gelation transition, 
Xp = 1/(1 — 9^) to macroscopic phase separation in the 
preparation ensemble and Xm = Xcrit to microphase sep- 
aration in the gel. 

In flg. [s] we show a phase diagram in the Xm-A* plane for 
the special case q = 0. [The spinodals of the asymmetric 
case can be recovered by replacing X by (1 — 9^)x.] The 
dashed line = I separates the gel state and the liquid 
state. The latter is further divided into a mixed and a 
macroscopically phase separated liquid at Xm = 1 (solid 
line). The dotted line Xm = Xcrit separates the mixed gel 
from the microphase separated one. 
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FIG. 3: Pseudo phase diagram of the polymer blend in 
the x-M'Pl^i^s- The state of the system is, however, history- 
dependent (see text for details). 



The diagram in fig. |3| is not a true equilibrium phase 
diagram, because the state of the system also depends 
on the preparation ensemble via Xp. In particular, the 
microphases are only obtained if crosslinking takes place 
in the homogeneously mixed phase. As a consequence, 
the transition line fi — 1 cannot be crossed from the 
macro- to the microphase separated state for Xm > 1- 
The history of the gel is indicated by a path in the dia- 
gram. Of particular interest are the three paths a, b and 
c. Path a amounts to crosslinking close to macroscopic 
phase separation, Xp — > 1, and the endpoint A corre- 
sponds to a homogeneously mixed gel having long-ranged 
frozen-in charge fluctuations. Along path b the system is 
crosslinked in a preparation state that is far away from 
macroscopic phase separation. The endpoint B corre- 
sponds to a homogeneously mixed and rather weak gel, 
just crosslinked enough to be solid-like. Along path c, the 
system is prepared in the same way as on path 6, how- 
ever more crosslinks are introduced, which strengthens 
the gel. Subsequently, the temperature is lowered (the 
incompatibility Xm increased) , so that the endpoint C is 
close to microphase separation. These three histories are 
representative in the following sense. Each endpoint cor- 
responds to a state close to one critical point, as discussed 
above: point A is close to macroscopic phase separation 
in the preparation ensemble, Xp = Ij point B is close to 
the gelation transition, /i = 1; and point C is close to mi- 
crophase separation, Xm = Xcrit- The three states of the 
system, corresponding to the endpoints, will be discussed 
in detail in the following sections. 



E. Charge density correlations in the mixed gel 

In this section we discuss the homogeneously mixed 
gel phase, for which the order parameter for phase sepa- 
ration vanishes. Nevertheless, there are thermal as well 
as quenched charge fluctuations on various length-scales. 



These can be detected with help of multiple correlation 
functions. On the Gaussian level of approximation, these 
correlation functions are given by the inverse of the Hes- 
sian matrix (p4|): 



n+l 



{A ^)aii 



(32) 



The correlator that is off-diagonal in replica space ac- 
counts for the frozen-in correlations and will be termed 
the glassy correlator. It is given by 



b{b+c) . 
c-Xi ' 



(33) 



see eq. ^ for the definition of ^. The replica-diagonal 
correlator is the scattering intensity 

and the variance (or connected correlator) is given by 

S^Uk) ■■= Ssc(k) ~ Sgiik) = ^ . (j- - Xm) • (35) 

Whereas the glassy correlator Sgi{k) describes the static, 
frozen-in correlations, the variance S'var(fc) quantifies the 
volatile, thermal fluctuations about the mean value. The 
scattering intensity 5sc(fc) is the sum of both contribu- 
tions and covers both, thermal and static charge inhomo- 
geneities. 

Restriction to symmetric blends 

In the following discussion of S'sc(fc), 'S'gi(fc) and Svarik) 
we shall confine ourselves to the case of symmetric blends, 
yet without loss of generality: The scattering functions of 
the asymmetric case are recovered via multiplication by 
7 = (l-g^) and the rescahngs Xp ^ 7Xp and Xm ^ 7Xm- 
Furthermore, the distance to phase separation is replaced 
with the distance to the spinodal in the asymmetric case. 
In the range between the equilibrium phase transition 
and the spinodal, the results then describe an under- 
cooled mixture. 

Length-scales 

The correlation functions are characterized by three 
length-scales which are determined by the parameters 
ip-i Xp: Xm) of preparation and measurement conditions: 

First, there is the typical localization length ^ of the 
monomers in the gel fraction, i.e. the mean mesh size of 
the gel. From eq. (^2|) we can infer that this length-scale 
is roughly given by 



6 := 1 / Vm^- 



(36) 
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Second, there is the decay length of the pre-critical 
demixing fluctuations prior to gelation. This approxi- 
mately reads 



1 / yi^- 



(37) 



The third length characterizes the pre-critical fluctu- 
ations of microphase separation, and is approximately 
given by 



Cm 1 / VXcrit(Ai) - Xr 



(38) 



The three length-scales measure, or are given by, the 
inverse distance to the phase transitions of gelation and 
demixing in the pre-crosslinking blend, and microphase 
separation in the gel; hence they grow large when ap- 
proaching their respective transitions. In the follow- 
ing, we shall essentially discuss three limiting regimes, 
in which the correlation functions are determined by one 
of the three length-scales 

• Cp > 6, Cm (point A), 

• 6 > Cp,Cm (point B), and 

• Cm > 6,Cp (point C). 



Glassy correlations 

The glassy correlation function Sgi{k) describes time- 
persistent charge inhomogeneities due to crosslinking. If 
the preparation ensemble is close to phase separation 
then instantaneous crosslinking will freeze in these fluctu- 
ations, and Sgi will be dominated by the pre-crosslinking 
fluctuations, giving rise to a high value at zero wave- 
vector. If, on the other hand, the preparation ensem- 
ble is in a well-mixed state then crosslinking will in- 
troduce completely random, static charge fluctuations, 
which subsequently can be enhanced by approaching the 
microphase separation transition in the gel. In the fol- 
lowing, we discuss the three limiting cases (i.e. points A, 
B, C) in detail. 

We first consider a gel that is prepared from a melt 
close to phase separation, i.e. ^ (6: Cm) ^ 1, corre- 
sponding to point A in fig. ^j. The network can freeze-in 
correlations on length-scales larger than or comparable to 
its mesh size. For ^p ^ ^i, the pre-crosslinking fluctua- 
tions have long enough scales to be frozen. Consequently, 
the glassy correlations reflect the pre-crosslinking fluctu- 
ations: 



c 



(39) 



The glassy correlations are proportional to ^p and de- 
cay on the scale A; ~ set by the fluctuations of the 
preparation ensemble. An example is included in fig. |4| 

In a weak gel, i.e. if ^ (CpjCm) ^ 1 (point B), the 
network is rather wide-meshed, so that the fluctuations 
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FIG. 4: Correlation functions for a gel prepared close to 
phase separation: Sgi, 5var and Sbc in units of 4^pCm/Cf for 
Cf = 10', = 10' andC,n = 10. 




FIG. 5: Correlation functions for a weak gel: Sgi, Svar and Ss, 
in units of for = 10^ = 10 and = 10. 



at preparation cannot be frozen in. Instead, there will 
be static charge fluctuations {^k) 7^ on the scale of the 
network, which are completely random and hence van- 
ish, if averaged over crosslink configurations, [(^'fc)] = 0. 
They do, however, contribute to the glassy correlations, 
which are given approximately by 



(40) 



These fluctuations always decay on the length-scale of 
localization, but they are enhanced in magnitude when 
approaching microphase separation. An example of the 
glassy correlations in this range is given in fig. |5|. 

The cross-over between the two scales is demonstrated 
in fig. |6[ which shows S'gi(A:)/S'gi(0) far from microphase 
separation for ^1 = 10. For the leftmost curve = 10^ ^ 
^2 = xoO, and hence the decay occurs at k ^ . Upon 
decreasing ^p, the curves shift to the right, until, for ^1 
^p, the decay is determined by ^1. The inset shows the 
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FIG. 6: Scale crossover in the glassy correlation function 5gi, 
normalized to the value at fc = 0: Crossover from preparation 
close to demixing to a weak gel, with Xm = 0.1, = 10^, 
and = 10^ 10''■^ . . . , 10^ (from left to right). Inset: Half- 
width fcjyj of Sgi{k), crossing over from to ^j"'^ (dashed 
line). 




FIG. 7: Correlation functions for a gel measured close to 
microphase separation: Sgi, Svai and Ssc in units of ^m/^f 
for = 10^, £,1 = 10 and = 10^. The wave-number 
squares are measured in units of k^. 



half- width at half- maximum as a function of f „ ^ . 

Close to microphase separation, i.e. ^ (6iCp) ^ 1 
(point C), Sgi{k) is dominated by the critical fluctua- 
tions towards microphase separation. In the absence 
of crosslinks there would be large-scale fluctuations to- 
wards macroscopic demixing. In the gel, displacements 
are bounded by the localization length, so that S'gi(A;) de- 
velops a peak at kc ~ ^, where A2(fc) becomes small; an 
example is included in fig. 0. Approaching the transition, 
the peak diverges as A^^(fc), and the glassy correlations 



can be approximated by 
Sgi{k) cx 



Q 



iXcM - Xm + - fc?)2w;"(l)/fc2)2 ' 

(41) 

where w" denotes the second derivative of the scaling 
function introduced below eq. ( p2|) and defined in ap- 
pendix 



Thermal fluctuations 

The variance S'var(fc) of the charge fluctuations is in- 
dependent of the conditions at the time of crosslinking. 
Hence, there are only two competing length-scales, and 
^m. In case of a weak gel, i.e. » ^ 1 (point B in 
fig. H), eq. ( |35| ) reduces to 



'5'var {k) 



Sn 



1 -I- k^il/i 



(42) 



for small fc, decaying with a half-width of ki/2 ~ "v/S/Cmj 
provided is not too small. Far away from the demix- 
ing threshold, the fluctuations are hardly affected by the 
network and look like critical fluctuations approaching 
macroscopic phase separation. An example is shown in 
fig. I 

Close to microphase separation, i.e. S> 3> 1 
(point C), the fluctuations grow with increasing k until k 
reaches the inverse localization length k ~ ^\^, beyond 
which they are strongly suppressed by the network. The 
variance is approximately given by 



'S'var (fc) 



A2(fc) 



Xcrit(M) - Xm + i(fc2 - klYw"{l)/kl 



(43) 



revealing a peak at fcn ^ 1 /^i that has a height propor- 
tional to f,^. See fig. M for an example. 



Scattering intensity 

The behavior of the scattering function in the vari- 
ous regimes can be inferred from the behaviors of Sgi{k) 
and S'var (fc), as Ssdk) is just the sum of them. A weak 
gel (point B) preserves only a small amount of the pre- 
crosslinking fluctuations and can hardly restrict thermal 
fluctuations. Hence, Ssc{k) ~ S'var (fc), so the scattering 
function decays on the scale k ~ ^nii see fig. |[ In the 
other regimes, thermal fluctuations are suppressed by the 
network, and Ssdk) w Sgi(fc). In a gel prepared close to 
phase separation (point A), the scattering function de- 
cays on the scale k ^ (fig. ^), whereas a gel measured 
close to microphase separation (path C) reveals a peak 
at fc ~ fco: diverging at the transition (fig. 
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V. MICROSTRUCTURES 

At 

Xia — Xcritj the homogeneous gel becomes unsta- 
ble with respect to phase separation. As we have seen 



in section IV C , the instability first occurs for nonzero 



wave-numbers, indicating that the gel undergoes micro- 
scopic rather than macroscopic phase separation. In this 
section we investigate various microstructures, such as 
hexagons and lamellae with a definite orientation, as well 
as a superposition of many random orientations. The se- 
lection of a particular microstructure depends sensitively 
on the compressibility and the charge imbalance. We 



first discuss the simplest case of incompressible, symmet- 
ric mixtures, and then go on to investigate the effects of 
charge imbalance and compressibility. 



Incompressible, symmetric mixtures 

Our analysis of microphase separation is based on the 
effective free energy ( |l5|) in the gel phase. We expand 
it around ^' = up to fourth order, in the presence of 
a nonzero gel order parameter Q given by eq. (|2^). The 
expansion reads 



ai2 k 



IE E' 



-|-fc2-|-fc3-|-fc4,0 



5D((fei+fc2)2) 3 



" ^1,2,3,4 

Here, A^^ a^i^) given in eqs. (24,2^), and the vertex of the fourth-order term is given by the Wiener correlators 



<73(fci,fc2)- / dsids2ds3 (e-'5:Li'=.'-(«.)^ 



w 



k3 = — ki—k2 



and 



g^4(fci,fc2,fc3) = j dsids2ds3ds4 (e~'^-=^^"'^^''''^j 



ki = -Y.l = ik^ 



r 



(45) 



(46) 



The fourth-order term in eq. apparently depends 
on the spatial structure of the microphases, and is re- 
sponsible for the pattern selection as well as for the wave- 
number selection in the microphase separated state. A 
particularly simple pattern are lamellae with sinusoidal 
modulations in real space: 



p 



0, for a = 0, 

V^{5p^k + '5p,-fe) V'; otherwise. 



(47) 



This ansatz is replica-symmetric, apart from the zeroth 
replica, which reflects the preparation ensemble. The 
wave-length 27: /k and the amplitude ij^ > Q arc vari- 
ational parameters subject to optimization. Insertion 
into (H) yields 



/({*}) ^ • lini .F„({*}) 



A2(fc2).^2^i54(fc')-^\ (48) 



with A2 defined in eq. ( P7| ) and 
1 f{gi{k,k)f 



= l-\e + 0{k^) (49) 



(see appendix g) . 

At the onset of the microphases, the amplitude goes to 
zero continuously and the optimal wave-number is given 
by fee [see eq. (|3l|)], implying a domain-size of the order 
of the localization length of the gel. Beyond the critical 
point, the amplitude is nonzero, and the wavenumber 
deviates from its critical value fee- Both are obtained 
from a variational optimization of the above free energy: 



V'min(fc) 



A2(fc) 



(50) 



and 



"Tnin ""c 



Co ■ (Xm - Xcrit) + 0((Xm " Xcrit)^) • (51) 



The constant cq can be computed analytically only for 
a weak gel, by further expansion in powers of /i — 1, 
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which yields Cq « 2.03(/x — 1). The optimal amplitude 
V'min(^min) grows continuously with Xm - Xcrit: 



2 

min 



Xr, 



Xcrit 



Xcrit 54(fc? 



+ O ((Xm - Xcrit)') 



(52) 



Other simple structures, such as a hexagonal stack of 
cylinders or a bcc crystal of spheres, have higher free 
energies. The same holds for a superposition of several 
sinusoidal modulations like eq. (p^), with the same wave- 
numbers but different directions. As we shall see below, 
these conclusions depend on the symmetry of the mixture 
and its incompressibility. 



Effects of asymmetry 

The most important effect of the asymmetry is a third- 
order term in the Landau expansion, rendering both the 



macrophase separation of the uncrosslinked liquid and 
the microphase separation of the gel first-order transi- 
tions. 



To keep the discussion simple, we neglect deviations 
of fc^ from and drop the A:-dependence of the higher- 
order terms in the Landau free energy. With the ansatz 
*g = (1 — Safi) corresponding to the mixed state in 
the preparation ensemble and replica-symmetric phase 
separation in the measurement ensemble, this leads to 
the free energy density 



l-3g^ 
12 



*fci*fe2*fe3*fc4 • '5fcl+fc2+fc3+fe4,0 , (53) 



with A2(fc) defined in eq. (P7|). 

Besides the lamellar microphases already discussed in 
section we now consider two additional morphologies: 
Cylindrical phases having parallel orientation, aligned 
on a honeycomb lattice in the perpendicular plane, and 
spherical domains on a body centered cubic lattice. Al- 
though a randomly crosslinked blend will probably re- 
veal only local order, the regular structures are useful for 
constructing a simple and tractable ansatz for the mi- 
crophase separated state: 



^ i—1 



(54) 



with m = 1 for lamella, m = 3 for hexagonally ordered 
cylinders, and m = 6 for spheres on a bcc lattice, and 
the corre spon ding lattice vectors {rii} being defined in 
appendix A 1 . 



With the lattice ansatz (p4 ), the evaluation of the 
higher-order sums in eq. (p3|) amounts to counting the 
number of possible "loops" of two, three and four lat- 
tice vectors that add to zero. This is carried out in ap- 
pendix Ala, yielding 'i'k'i'-k = 2^?^, independent 



of the morphology, and also 

ffei^fcs^fca ■ 4i+fc2+fc3,0 = 4'"^*^ (55) 



fel,2,3 



Vl/fc^*fc2*fe3*fc4-4i+fc2+fc3+fc4,0=cf^*', (56) 



fcl,2,3,4 

where 

'-a 



= 0, 
= 6, 



„(3) 



= 10, 



ci^) = 15. 



(57) 
(58) 



Thus, the free-energy density becomes 



l-g2 



+ 2q^ 



2 , (l-3g')4"^' 



12 



(59) 



Here, A2 has been evaluated at because (in this section) 
we are not considering deviations of the wave-number 
from its critical value. 

For TO = 1, the third-order term vanishes, even iiq ^ Q, 
so the transition remains second order and the spinodal 
indeed indicates the equilibrium phase transition point 
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with respect to lamellae. In contrast, for cylinders and 
bcc spheres the equilibrium transition point Xt is shifted 
according to 



1 



(1 - q^)Xt (1 - q^)Xc 



(9crV3)' 



4(2^2 + (1 _ 3g2)4"Vl2) 

/ (m)\2 

72q^ + 3(1 - 3g2)4'") 

8q^ _ 24 2 4^ for m 3 

45 - 27g2 ~ 135 ^ ^ ^'^ ' (cylinders) 

32g^ _ 32 2 g( 4^ form = 6 
135 - 189g2 135 ^ ^ (bcc spheres) 

(60) 

In the asymmetric case, the bcc spheres yield the lowest 
equilibrium transition point of the three possibilities con- 
sidered, i.e. microphases first occur with bcc symmetry. 
This is to be expected, as the ratio of surface to volume of 
the minority phases is minimal for spheres embedded in 
the majority phase, and is in agreement with the finding 
of Alexander and McTague |]13| of a general preference 
for bcc symmetry in crystal nucleation. Note, however, 
that the Landau expansion is only valid for small q^, for 
which the transition is weakly first order. In particular, 
for g2 — 5/7^ where the right hand side in the last line of 
eq. ( |60|) diverges, the Landau expansion breaks down. 



Effects of compressibility 

A compressible system can avoid unfavorable A-B con- 
tacts and lower its energy by diluting mixed regions hav- 
ing many such contacts, which are characterized by a 
small absolute charge density, and condensing regions 
that are rich in either A or B, which have a high absolute 
charge density. Mathematically, this becomes apparent 
via a nonzero value of the saddle-point of the density p 
in a phase separated state. For simplicity, we restrict 
the discussion of compressibility effects to the symmet- 
ric case, where the shifted and the original density fields 
coincide. In this case, the saddle point of p is given by 



Pk 



2(l/A"+5D(fc2)) 



E 

fcl,2 



n,n.-Sk+k,+k,,o. (61) 



For the simple example of lamellar microphases described 
by a single wave- vector k as in eq. (^ , eq. ( |6l]) predicts 
density-field modulations having wave-vector fci = zt2k, 
i.e. with twice the wave- number of the charge-density 
modulations. This is intuitively clear: Along one spatial 
period of the charge-density modulations, their modulus 
or square, and thus the mass density, oscillates twice. 




FIG. 8: Coupling of mass (p) and charge (^P) density. Zones 
of large charge-density modulation are condensed, hence the 
mass-density is modulated with half the wave-length of the 
charge-density modulations. 



corresponding to the half wave-length or the double wave- 
number; this is illustrated in fig. g. 

The compressibility is controlled by the strength of 
the excluded volume interaction, i.e. the parameter 
A = Am — p/V^. To study microphase separation in 
the symmetric but compressible case, we integrate out 
the density field, keeping A finite. With the ansatz 
= (1 - 5a,o)*, we obtain 



+ 



\ k 
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2 Aoff , 



fel,2,3.4 



X Ski+k2+k3+k4,0 ■ (62) 



Here, Xes = X ~ p + 1, and we have dropped the k de- 
pendence in the higher-order vertices, thereby restricting 
the domain size to its critical value, determined by fcc or 
the localization length of the gel. 

To account for the potential randomness of the 
microphase pattern, we extend the previous lamellar 
ansatz by allowing a superposition of Z one-dimensional 
waves, each with identical wave-number but random 
phases $2 and wave- vector orientations n^, i.e., 



~Z 



z 

E 



(63) 



corresponding to 2"^ / [V Z ' ) -X^z^i *^os {kcUzX + $2) in 
real space. The optimal number of orientations will be 
determined later. A few examples of such random mor- 
phologies are shown in fig. ^ the number of phases rang- 
ing from 1 to 100. 

We assume that none of the orientations are 
coUinear, and thus the quadratic sums in eq. (|62|) yield 



E 

k 



A2(fcc)*fc1'-fc = 2A2(fce)4'2 
^'«'fc*-fc = 2*2. 



and 



(64) 
(65) 
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FIG. 9: Superposition of lamellae having random orientations 
in real space in two dimensions. The pictures show an area of 
10 by 10 wave-lengths, and the local amplitude is indicated 
by the grey scale in arbitrary units for Z = 1,5, 10 and 100 
orientations. 

To compute the fourth-order sum in eq. ( |6^ ) , we have to 
count the number of possible closed loops of orientations. 
Because of the randomness, the existence of quadruples 
of orientations able to form a closed loop is very unlikely, 
except for the degenerate planar case of pairs of opposite 
vectors (iriz, ztn^/), and hence we disregard non- planar 
loops. Single orientations allow for the construction of 
quadruples (n^, n^, — n^, — 712) that can be ordered in 
(2) = 6 ways. Quadruples (in^ , in^' ) of two pairs of 
different orientations can be ordered in 4! = 24 different 
ways, and there ar e \ Z{Z — 1) such pairs. Thus, the 
quartic sum in eq. ( p^ ) yields 

*fci*fc2*fc3*fc4 • '^fcl-|-fc2-|-fc3+fc4,0 

'=1,2,3,4 

YlZiZ -\)+%Z 4 / 1 \ 4 , ^ 

— ^ -^^\^-2zr - (66) 

Inserting the sums into the free-energy density we ob- 
tain 

The fourth-order term depends on Aeff and the number 
of components Z . It has to be positive to guarantee sta- 
bility, and this requires Aoff > 1. The sign of the 0{Z~^^ 
term determines the optimal number of random orienta- 
tions. For low compressibilities, i.e. for AeS > 3/2, the 
term is negative and the free energy grows with an in- 
creasing number of orientations, hence the simple lamel- 
lar morphology is favored. For a rather compressible sys- 
tem having, in contrast, 1 < Acg < 3/2, the effective free 
energy decreases with increasing Z, favoring an "infinite" 
number of orientations and hence a random pattern. 

VI. CONCLUSIONS 

In this paper we have analyzed a microscopic model of 
crosslinked polymer blends, built on the Edwards model 
for a polymer melt and generalized to two components, 
which are mutually incompatible. Random crosslinks are 
introduced according to the Deam-Edwards distribution. 



also generalized to include concentration fluctuations at 
the instant of crosslinking. Thereby, the concentration 
fluctuations in the melt are partially frozen in and sus- 
tained in the gel phase. Apart from these correlations, 
which are present at preparation, the crosslinks are taken 
to connect monomers irrespective of their charge. Hence, 
within mean-field theory, the resulting gel is identical to 
the one made from just one species of polymer. However, 
concentration fluctuations are present, and have been 
computed on the Gaussian level of approximation. Of 
particular interest are the frozen-in or glassy fluctuations, 
which reflect the preparation state. In general, the net- 
work can only quench fluctuations on length-scales larger 
than its own mesh size, which is roughly given by the lo- 
calization length of mean-fleld theory. If the preparation 
ensemble is close to macroscopic phase separation then 
the length-scale of these fluctuations is large compared to 
the mesh size, so that the glassy fluctuations are given by 
the concentration fluctuations in the preparation state. 
If, on the other hand, the preparation state is far from 
phase separation, the frozen-in charge fluctuations follow 
the network pattern, and hence are completely random, 
because the crosslinks are not sensitive to the species. 
The thermal concentration fluctuations are independent 
of the preparation state. 

Lowering the temperature in the gel, or equivalently in- 
creasing the incompatibility of the two species, gives rise 
to phase separation, the spatial extent of which is lim- 
ited by the mesh size of the network. The length-scale of 
the resulting "micro"-phases can thus range from almost 
microscopic to nearly macroscopic scales, depending on 
the degree of crosslinking of the gel. The instability to- 
wards microphase separation is signaled by a divergence 
of the time-persistent as well as the thermal fluctuations. 
The emergent microstructure is shown to depend sensi- 
tively on charge imbalance and compressibility. The lat- 
ter allows for random patterns with a unique wave-length, 
i.e., the localization length, whereas in the incompress- 
ible system lamellae are favorable for balanced mixtures 
and hexagonal patterns for imbalanced mixtures. 

We now compare our results to previous phenomcno- 
logical approaches, many of which have focused on the 
issue of microphase separation, de Gennes argued that 
the charges in the crosslinked gel cannot move freely but 
are displaced in analogy to the charges in a dielectric 
material. He introduced a polarization P which, as in 
electrostatics, is determined by the charges according to 
V • P = — 5*. In the limit of weak segregation, the free en- 
ergy is quadratic in the polarization, and is simply added 
to the free energy of charge fluctuations, resulting in 

/ ({*}) = ^ E' (^-t - X- + + ^^^-^ ' 

(68) 

where C is a coefficient of "internal rigidity". The 
above free energy leads to an instability at finite 
wave- number (microphase separation), but predicts that 
limfc^o 'S'sc(fc) = 0, in disagreement with experiment. 
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The nonzero scattering intensity at zero wave-number is 
due to the frozen-in charge fluctuations present at prepa- 
ration. To account for these fluctuations, Benhamou et 
al. 1^ refined the analogy to a dielectric by including a 
Debye-Hiickel screening of the "charges", which permits 
long-range inhomogeneities leading to a non-vanishing 
zero-angle scattering. The screening length k is deter- 
mined self-consistently, by assuming that the scattering 
intensity at A: = is not affected by the crosslinking as 
long as the temperature remains unchanged after prepa- 



ration |14| . The free energy in the quadratic approxima- 
tion then reads: 

/ ({*}) = I E' {Xcr- - X.n + e + . 

(69) 

This expression can be compared with eq. ( p3[ ) in the 
quadratic approximation: 



/ ({*}) = ^ E' f — - 9D{k') + m(m - l)^(fc'(M - 1))) ^k^-k 

^ ^ E' (i - + + - - 1))) 



(70) 
(71) 



In the last line, we have expanded the Debye function for 
small wave-number. We see that the microscopic model 
indeed agrees with phenomenological theories, provided 
we identify the phenomenological terms with the order 
parameter of the gel. The wave-number dependence of 
the order parameter is not a Lorentzian; nevertheless, 
it decays monotonically with k, the relevant length-scale 
being given by the localization length. Hence, the some- 
what mysterious screening length is unambiguously iden- 
tified with the localization length, which is computed self- 
consistently. Thereby, the microscopic model substanti- 
ates the picture of de Gennes and, furthermore, allows 
the computation of the parameters and functions that 
are beyond the phenomenological approach. 

The frozen-in fluctuations were first addressed by Read 
et al. 0, who considered a blend of polymer chains an- 
chored at both ends to randomly chosen fixed points in 
space, in order to account approximately for the local- 
ization of chains due to the crosslinks. Read et al. make 
reasonable but ad hoc assumptions about the distribution 
of the quenched random end-to-end vectors, and solve 
the resulting model within the random phase approxi- 
mation. They obtain a scattering function that exhibits 
a nonzero value at fc = 0, due to the random, quenched 
fluctuations. In addition they compute the thermal, as 
well as the glassy, charge fluctuations: 



limit of small wave-number 

"5'var (fc) 



i-xm + k^^ + xm(M - iMkyifi - 1)) 



(1 - Xm + fc2 2^ + Xm(A^ - \)w(k-^l{il - \)W 



Both approaches, the phenomenological one and the mi- 
croscopic model, predict a divergence as microphase sep- 
aration is approached, with the glassy correlations di- 
verging twice as strongly as the thermal ones. 

The work reported in the present paper can be ex- 
tended in several directions. First, we have worked only 
on the level of mean-field and Gaussian fiuctuations. It is 
known that the microphase separation transition in the 
symmetric case is rendered first order by fluctuations , 
and hence it would be interesting to see the effect of fluc- 
tuations, even though the critical region is expected to 
be small [|l6|. Another extension is a crosslink proba- 
bility that depends on the species. This would allow us 
to study, among other things, interpenetrating networks. 
Finally, it would be interesting to look at the dynamics 
of microphase separation. 



APPENDIX A: MICROPHASE MORPHOLOGY 



'S'var(fc) — 



1 



5g,(fc) 



Xcrit - X + + ^ 



(Xc 



X 



C ^2 



and 



(72) 
(73) 



where po{k) is the frozen-in concentration. The above 
results are in close correspondence to the results of our 
analysis, presented in eqs. (p3|,p5|) and evaluated in the 



To investigate microphase transition, we assume that 
the phase-separation pattern can be described by a first- 
harmonic ansatz having a dominant wave-number k and 
a definite lattice structure: 



*fe' = ^' y— 7 ^ {Sk',+kni + 5k' -km) , 



(Al) 



with lattice vectors Ui 
amplitude 5*. 



{rii \ i = 1 . . . m} and an 
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Lattice structures 



We consider three particular morphologies, which are 
known to occur in the microphase separation of regular 
copolymer melts 18 : 



lamellae {m — 1): Alternating sheets rich in A and B; 
onc-dimensional order. Lattice vector: rii — 

(1,0, or. 

cylinders {m — 3): Close- packed, i.e. hexagonally ar- 
ranged, cylindrical domains, A in B or vice versa; 
two-dimensional order. Lattice vectors: rii — 
(1,0,0)^, n2 = (-l/2,V3/2,0)^ and ng - 
(-l/2,-V3/2,0r. 

bee spheres (m = 6): Spherical A-rich domains in B, 
or vice versa, on a bcc lattice in real space; 
three-dimensional order. Lattice vectors of 
the corresponding fee lattice in Fourier space: 



m = (i,i,or/V2, 
(1,0,1)^/72, m = 

(-l,l,0)'^/\/2 and Tie 



712 = (0,l,ir/%/2, 
= (1,0,-1)VV2 ns 
(0,-l,l^/^/2. 



Note that a e G a ^ G; therefore we introduce 
the symmetrized set of lattice vectors, (G+ := {n | n G 
G V n E G}. For to > 1, the set G of lattice vectors is 
not minimal in the sense of linear independence: for any 
two vectors a ^ b E G, the difference a — b is included 
in G+. Rather, the vectors are chosen such that the ki E 
G^ point towards the directions of all nearest-neighbor 
lattice sites. 



a. Wave-vector sums 



Inserting the ansatz (Al) into the expansion of the 
Landau free energy of the random copolymer melt or the 
crosslinked homopolymer blend yields sums of the type 



Y,' /2(fei, . . . , kp-i)^k, ■ ■ ■ 'fkJY.l^i 



W -VP A / \ 



VP 



(2m)P/2 



^ /2(fei,...,fcp_i)4i+...+fcp,o- (A2) 



fci „gg+ 



In the quadratic sum the vertex function can be fac- 
tored out, so that 



2to 



loop type 



lamellae cylinders bcc spheres 



2- loop 2 6 12 

3- loop 12 48 

4- loop planar 6 90 396 
nonplanar - - 144 
total 6 90 540 

TABLE I: Number of closed loops of p lattice vectors ("p- 



^ ^ <Jfci+fc2,o ■ f{k^) ■ (A^3) loops") for different morphologies. The 4-loops are divided 



fci,2eG+ 



into intra- and extra-planar loops. 



The higher-order sums, however, strongly depend on the 
morphology of the microphases. In the simplest case, 
/(fci, . . .) = 1, the computation of these loops amounts 
to counting the number of closed loops that can be con- 
structed with the vectors in G"*". In general, the loops 
must also be classified with respect to their shape, i.e. 
planar or non-planar, as distinct shapes yield distinct 
values of the vertex functions. The counting and clas- 
sification have been done, e.g., in p8| , with the results 
shown in table Ala. For /(fc, . . .) — 1, the third- and 



fourth-order sums are given by 



1 



for TO = 1 (lamellae), 

• <^ for TO = 3 (cyhnders), (A4) 

2/V3 for m = 6 (bcc spheres). 
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and 



3/2 for TO = 1 (lamellae), 
• <( 5/2 for TO = 3 (cylinders), (A5) 
15/4 for TO = 6 (bcc spheres). 



APPENDIX B: VERTEX FUNCTIONS 

In the following, we define and compute the vertex 
functions appearing in the Landau expansion of the ef- 
fective free energy, which are integrals over Wiener cor- 
relators of the type 



w 



^ C=i ^ 



cc 



a derivation of cq. (Bl) can be found in (l|. The cor- 
relator vanishes unless the wave-vectors sum to zero in 
each replica. If just single- replica quantities are involved, 
the correlator factorizes. 



for pairwise distinct ai, . . . , ttm, where 
the unreplicated Wiener average. 



\W 



denotes 



1. Definition of the vertex functions 



The second-order coefficients of the Landau expansion 
are governed by the Debye function 



5d( 



(fc2) := / dsids2(e-"=('-(^^)-'-(^^») 



1 , . >w 

(si)-f (S2)) \ 

/ n 



= I dsids2 (e 



the scattering function for a non-interacting Gaussian 
chain. 



The third-order correlators, 



exp<^ -i XI '^Ci • r^HsCi) - ■■■ 

-i J2 K.-r'^-i^cJ 

= /expj - i X hir{sci) 
^ ^ Ci=i 



Cn. = l 
V W 



w 



exp <j - i X ^C".'^(*C™) 



w 



S'3(fcl,fc2) 

:= dsids2ds3 {^e-' 2 



w 



dsids2ds3 (e-'S'=i'^-^("-) 



w 



/C3 = — fel— fc2 



and 



5*25^(^1,^2):= / dsids2ds3(e 



ifci(r(si)-r(s3)) \ / -ifc2(i'(s2)-''(s3)) \ 



W 



= y dsids2ds3 ^e-"=i''°'(^i)-*=^''°'("^)-"=^('*^) 



V 



describe the correlation between three one-replica fields 
and the correlation between two one-replica-fields and a 
higher-order replica field, respectively. 



Finally, the fourth-order correlator is given by 



dsids2ds3ds4 ^e^'S-i 



w 
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a. Lamellar case quire the "lamellar" case, in which all wave-vectors are 

coUinear. We note that 53(fc, —k) — gY){k^), and define 

The third- and fourth-order correlators depend on the 
directions of the wave- vectors. In particular, we re- 

I 



and 



,2, X -(e" - 1 + 4fc^ - 8/c'') 64(e-''' - 1 + fc^ - ifc*) 2(0"'' - 1 + P) 



9z{k ) := gaik,k) = 



12/c6 /c4 



/,2^ ,x -C-2fe- + (8 + 2P)c-fc + 4P ~ 7 , 2,2 17,4 .^/,6n 



,,2x /, , ,x 144fc4 -60fc2(9 + 4c-M + 784(1 - c-M - (l-c-4fe 



IS/fcS 



APPENDIX C: SCALING FUNCTION FOR THE 
GELATION ORDER PARAMETER 

The localization lengths t of monomers in the gel frac- 
tion of a crosslinked homopolymer melt or blend are 
distr ibute d according to the distribution p(t); see sec- 
The fraction of the gel and the distribution 



IV B 



tion 

arc determined from the self-consistent solution of the 
saddle-point equations with the order parameter hypoth- 
esis (pi|). The gelation order parameter, essentially the 
Laplace transform of p{t), is proportional to a scaling 
function uj{x), which is computed in fl^ in the asymp- 
totic regimes of small and large x. For convenience, we 
define a rescaled version of uj, 



'j{k^) := 2-w(V8fc2/3), 



(CI) 



This also absorbs a factor of two arising from the different 
length-scale used in (the Wiener Hamiltonian used 
therein differs by a factor of two). 



1. Interpolation formula 

The scaling function uj{k) defined in can be de- 
scribed asymptotically by 

f uji{x) := 1 - O.U09x^ + O.lSlGx"^ for x < 1, 



^8- 1.678^ 

x(l + 
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40\/2-1.678; 



for X > 1. 



In order to access the whole range of < x < oo we 
interpolate between the asymptotic regimes using the in- 



1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



2 
X 



FIG. 10: Scaling function ij{k) versus k: Asymptotic expres- 
sions and interpolated function. 



terpolation formula 



Lu{k) 



' UJi{x), for X < ^, 
uJip{x), for i < a; < 2, 
uj2{x), for X > 2, 



(C2) 



with the interpolating rational function 

bo + bix 



1 + aix + a2X^ + a^x^ 



(C3) 



The coefficients ai = —0.055, 02 = 0.165, 03 = 0.139, 
60 = 1.023 and 61 = —0.194 are chosen such that the 
value and first derivative of (^ip(x) coincide with those of 
ll}i{x) at X = 1/2 and with those of 0^2(3;) at x = 2; an 
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additional sampling point is the numerical value u){x = 
1) = 0.664. 
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